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Abstract 

In this paper we discuss general tridiagonal matrix models which are natural extensions of the ones given 
in [3] and [I] . We prove here the convergence of the distribution of the eigenvalues and compute the limiting 
distributions in some particular cases. We also discuss the limit of fluctuations, which, in a general context, turn 
out to be Gaussian. For the case of several random matrices, we prove the convergence of the joint moments 
and the convergence of the fluctuations to a Gaussian family. 

The methods involved are based on an elementary result on sequences of real numbers and a judicious counting 
of levels of paths. 

1 Introduction 



Tridiagonalization is a standard procedure in numerical analysis. The advantage of tridiagonalization is that the 
eigenvalues do not change under this procedure on one hand and on the other hand the tridiagonal matrix is 
easier to study, both numerically and theoretically. 

The well known GOE, GUE and GSE random matrix models (see [7] for a standard reference), have the 
eigenvalue distribution given by the density 
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for (3 — 1,2,4 and Z n a is the corresponding normalization constant. 

For j3 = 2, tridiagonalizing the GUE ensembles, in [5] and [I], the authors arrive at 
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where all entries are independent and Xr is the x distribution with r degrees of freedom. Since the tridiago- 
nalization does not change the eigenvalue distribution, it follows that for this model the eigenvalues have the 
distribution given by (jl.ip . Moreover it turns out that for any arbitrary fi > 0, the eigenvalue distribution of 
the model (TO)) is given by (jl.ip . 

Obviously the models (|1.2[) are less complex and consequently one should be able to take advantage of this, 
particularly in the case of computations of expectations of traces of powers. In [3] and [4] the limit distribution 
and the fluctuations are studied. However some of the arguments used there rely on the particular form of the 
model and it's not clear weather these particular properties are really needed for the convergence and fluctuations. 
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Another model which is discussed in the literature is the Wigner ensemble which appeared for the first time 
in [17) and [16] . These are symmetric random matrices with upper diagonal entries independent of one another 
with mean zero and the same variance. For these ensembles, Wigner himself proved a form of convergence of 
the distribution of eigenvalues to the semicircle law. The main method available here to study the limiting 
eigenvalue distribution and fluctuations is so called moment method which consists in expanding the traces of 
powers and counting the contributing terms. There are various sources using this method, among the so many 
we mention for instance [10] and the survey paper [5] for various combinatorial but also analytic approaches. 
For the problem of fluctuations from the limiting distribution, a very general form can be found in pQ . Another 
use of the moment problem is in [9] for universality at the edge of the spectrum. 

In the context of tridiagonal models we would like to introduce and discuss the analog of the Wigner ensembles 
and prove the convergence of the distribution of the eigenvalues and the fluctuations using the method of 
moments. We show a nice and clean combinatorial way of doing this. 

At first, these models may seem to be an extension in form only. There are many reason we want to study 
these. The first one is that these seem to be the natural analog of the Wigner ensembles for the tridiagonal 
ensembles. It turns out that these ensembles obey nice properties as convergence of the empirical distribution 
of the eigenvalues and the fluctuations converging to a Gaussian family. Thus these can be seen as another 
universality property. The second reason, the main one is that tridiagonalization of a Wiegner ensemble outputs 
a random tridiagonal matrix. We are still far from understanding these tridiagonal matrices due to the fact that 
the entries of the resulting matrix are no longer independent. What we try here is to study models in which the 
entries are independent with the hope that these will shed light on the more intricate case with dependencies. 
The third reason is connected to the following problem. Take a band matrix of width growing with the size of 
the matrix. These models have been studied in the literature in some situations, but there are various cases 
where not much is known. Such case is the one in which the band width is the square root of the matrix size. 
To the knowledge of the author it appears that the convergence of the empirical distribution is not known. The 
tridiagonalization of such an ensemble produces a matrix whose entries have strong dependencies but we believe 
that studying these will bring to light some interesting phenomena. 

Our main matrix model in this paper is given by 
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where the entries are independent random variables. In particular if {d n }™ =1 is a sequence of iid normal random 
variables and b n — Xnp/Vfi, then we get (ll.2[) . One of the main properties used in [3] and [I] to study the 
limiting eigenvalue distribution and the fluctuation is the simple fact that Xr ~ \fr converges in distribution to 
iV(0, 1/2). Rephrased, it implies that in distribution sense 

lim b n /y/n = I- (1-4) 

n — >oo 

This together with the fact that d n are iid with finite moments, turn out to be sufficient for proving the 
convergence of the eigenvalues to the semicircle law for the rescaled matrix X n = -^A n . 

In what follows, for any matrix Y = [yi,j)2j—i, we use Tr n (Y) ~ Y^i=i Vi,i f° r the full trace and tv n (Y) = 
— Vi,i f° r the reduced trace. 

To outline the idea of this paper in one instance, namely the convergence in moments of the eigenvalue 
distribution, let's take the trace of the fourth moment of X n , which is 

tr n (A n ) — — ^ ] Oil ,*2 °»2 ,*3 °»3 »*4 °«4 .*1 • 

l<2l ,22 ,23 ,14 <n 

We want to show that this converges. Here Ojj are the entries of the matrix A n . Now since the matrix A n 
is tridiagonal, these terms are zero for \i u — i u +\\ > 2 for 1 < u < 4 with i§ = i\. Hence the only nonzero 
contribution is given by the sequences 12, 23, 14) with \i u — i u +i\ < 1- Let's call these sequences admissible. 
Now we rewrite 
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Since the indices ?i,?2,*3,*4 are in within finite distance from one another, for p larger than 3, the sum 
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depends only on p and not on n and 
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In the limit, for large n, one can ignore from the above sum, the terms E[Si], E[S < 2] and E[Ss] or any finite 
number of them. The key to our computations is the following simple result on sequences: 
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Applying this to = E[5 P ], one reduces the computation of the limit of E[tr„(X^)] to 

lim ^E[Sy. 

p — >oo p z 

Next, we notice that i\ = ji+p, ii = 32 +P, h = h +P and U = H+V, with \ j u -j u+ i\ < 1, 1 < U < 4, j 5 = j\ 
and rewrite 
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and the limit of E[5 p ]/p 2 reduces to the ones of the form 



lim Efajj+pj^+pa^+pjs+pCtjs+p.^+pa^+jjj'i+p] 



(*) 



If j« = ju+i, then aj„+p,i„ + i+p = dj„ +p , while for the case |j„ - j u +i \ = 1, we have a,j u+p j u+1+p = hence 
if at least one of the situations j u = j u +i occur, (|1.4[) forces the limit in (*) to be 0. This means that the only 
contributing terms are those with all consecutive j'a different. This happens if and only if (ji, j 2l ja, ji, ji) is one 
of the following 



71 = (0,-1,0,-1,0), 
74 = (-1,-2,-1,0,-1) 



72 = (0,-1,-2,-1,0), 
75 = (-1,0,-1,-2,-1), 



73 = (-1,0,-1,0,-1), 
76 = (-2,-1,0,-1,-2). 



(1.6) 
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For each of these strings, one can compute the limit. For instance, in the case of 71 = (0,-1,0,-1,0), 
according to l|1.4jl . 

lim -^E[aj 1+ p j 2+p aj 2+ p j 3+ paj 3+ p j i+p aj i+ p j 1+ p] = lim — E[6p_ x ] = 1. 

p — >oo p^ p — >oo p^ 

Similarly we get 1 for all the other terms corresponding to these 6 strings, therefore one gets 

lim E[ix n (X*)] = 2, 



which is the fourth moment of the semicircular law ^1[_2,2] (^)v4 — x 2 dx. 



For fluctuations, the general statement is give in Theorem [31 To be in tune with the convergence discussed 
thus far, we want to show how one can deal with fluctuations on the following calculation: 

lim E[(Tr„(x4)-E[Tr„(^)) 2 ] . 

n — >oo 

Using 11.51 we can write (neglecting a finite number of terms) 
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Combining the independence of the entries with the expectation ,we obtain that S v , q = if \p — q\ > 2. Notice 
also that S p ^ q depends only on p and q for \p — q\ < 1 and not on n. Therefore, using the following elementary 
fact 

lim — ^ = M =>• lim — - > a; p = — , 

p^oo p-^ n— >oo ri — ' 4 

it suffices to deal with 



p — >oo p- 

P=PO 



lim % 

p — >oo p^ 

for g = porg = p±l. We can assume that p > q, otherwise we can simply reverse the order of p and q for the 
following argument. Using the same route as before for the computation of the limit, it suffices to find 

ml1 ~3 (^[ a jl+Pj2+P a j2+Pj3+P a j3+P,j4+P a j4+PJl+P a j;+PJ 2 +P^ 
P p 

— E[aj 1+ pj 2+ paj 2+ pj 3+p aj 3+ pj 4+ paj 4 ^ (1-7) 

where the strings here Ji, J2, J3,j4 and j{, j' 2 , j'4 are admissible with max(j u ,j' v ,u = 1 . . .4, v = 1 . . . 4) = 0. 

Any appearance of equal consecutive indices in the sequence Ji, J2, J3, J4, ji or j[, f 2 , j' 3 , j'4, j[, forces another 
appearance of consecutive indices and such an occurrence means that at least two of the entries in the above 
limit contains diagonal term. Since the diagonal terms are bounded (in moments) , this forces the above limit to 
be 0. 

This implies that the entries in the limit (|1.7|) are from the subdiagonals only and (ji, J2, J3, J4, Ji) is one 
of the paths -f u in (| 1 . 6|) and j' 2 , j' 3 , f 4 , j[) is one of the paths 7 U or j u — 1. Let's assume for simplicity that 
we deal with the model in which b n has a Xnff/Vft distribution. Using the fact that \ r — ^fr ~ N(0, 1/2) in 
distribution and moments sense, one gets that for fixed k,l > 1, cov(6p,6p) ~ klp( k+l ~ 2 ^ 2 /3~( k+l ^ 2 /4. Using 
this, one can compute the limit above. For example, if (ji, J2, J3, J4, ji) = 72 and f 2 , j'4, j[) =73 — 1, then 
the limit in (jTTTJ) is 



lim 1 (E^JE^] _ E^jE^jE^]) = l im Ie^cov^^, b 2 p _ x ) = 2/T 3 . 

P *oo p p *oc p 



Using this argument combined with the judicious counting of the paths, one can prove that the fluctuations 
converge to a Gaussian family. 

Let us return now to the convergence of the empirical distribution of the model 11.31 Since there is nothing 
sacrosanct about (|1.4p . we can replace it by 

lim E[(6„/Vn) fc ] = m k , 

n — >oo 

where is a given number. Loosely speaking this says that b n /\/n converges in distribution to a random 
variable Y with moments given by m^. In this case one gets for 71 = (0, —1, 0, — 1, 0) that 

lim — E[aj 1+p j 2+ paj 2+ pj 3+p aj 3+ p j 4+p a,j 4+p j 1+p ] = lim —E[bp_ L ] = rn 4 

p — >oo p^ p — >oo p^ 1 
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and in general, collecting all terms, one gets 
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The contribution of the paths is as follows: 
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For example 71 crosses the line —1/2 exactly 4 times and that corresponds to the index 4 in 7714, while the path 
72 crosses the lines —1/2 and —3/2 twice, each of these giving an 7712 term with the total contribution being the 
product of these, namely mf. 

Here we note that the scaling i/n in (|1.4|) is not essential for the argument. A more general treatment is one 
in which y/n is replaced by n a with a > 0, and on this line of ideas the first result we prove is Theorem [T] in 
Section [2] which concerns the convergence of the traces of powers, both in expectation and almost surely. There 
is also a combinatorial relationship between the moments of the limiting distribution and the moments of the 
limit b n /n a via counting the number of level crossing for paths. We specialize the limiting distribution in the 
case b n /n a coverges to 1. As opposed to the Wigner ensembles we get here different distributions depending 
on the scaling used and in some cases even an explicit formula. In Proposition [1] we give examples of limiting 
distributions for the case b n /n a converges to a Bernoulli random variable. Also worth mentioning here is the 
fact that the limiting distribution can be described as the distribution in a certain sense of a random Jacobi 
operator. At the end of Section [5] we also discuss the first order deviation of the expectation of the moments of 
the distribution of eigenvalues. 

The convergence of the fluctuations is fully discussed in Section [3J Under the appropriate conditions and 
after properly scaled, the family {tr n (X^) — E[tr„(X„)]}fe>i is shown to converge to a Gaussian family where 
the covariance can be computed. 

In Section[4]we extend Theorems[T]and[3]to the cases of multiple tridiagonal random matrices. This resembles 
very much the framework of free probability distribution (see |15j for basics and more) and also the second order 
freeness discovered by Speicher and Mingo in |8J . The interesting part would be to define some kind of cumulant 
similar to the classical cumulants or to the free cumulant (cf. [TT]) and then define some sort of "independence" 
via properties of cumulants. 

Section [5] gives various situations in which the same arguments can be employed to extend Theorems [T] and 
[3] As a particular case is the band diagonal and an eventual extension to the case in which the entries of the 
matrix are not independent. 

The combinatorics in this paper is one for Dyke paths. In the case of Wigner ensembles there is another 
combinatorial approach by counting planar graphs as it is done in [T] . It would be interesting to see the connection 
between these two combinatorial methods, though we do not have a clear way of bringing them together. 

The study of tridiagonal models in which dependence of the entries is allowed is very important. This is 
motivated for once by studying Wigner ensembles via tridiagonalization. In this case the independence of the 
entries of the tridiagonal model is in general lost. We hope that further study of the tridiagonal models may 
turn useful in the study of other random matrix models as for example band models where the band width grows 
with the dimension. 

The general belief is that the tridiagonal random matrices are easier to understand. This paper is a materi- 
alization of this belief in one instance, the case when the entries are independent. We hope that further study 
will turn this belief into a scientific fact. 

Acknowledgements Special thanks go to Elton Hsu for his interest in this problem and valuable discussions 
during the preparation of this paper. Also many thanks to the referees for their valuable comments which led 
to improvements of this paper. 

2 Convergence of the Distribution of Eigenvalues 

Our approach is combinatorial and as such we will deal with the convergence of the distribution of eigenvalues 
from the moment points of view. Thus, for a matrix A, the fcth moment of the empirical distribution of eigenvalues 
is given by the trace of A k . Hence, we reduce the study of the convergence of the moments of the distribution 
of eigenvalues to the convergence problem of traces of powers of the matrix. 

We start with the following elementary lemma which will be repeatedly used in this paper. 
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Lemma 1. 1. If x n is a sequence of real numbers, then for any p a > 1, and s > — 1, we have 

Hmsup M < M =► limsup * £ \**\ * TZ\' ^ 

P=Po 

I?. I/x„ is a sequence of real numbers, then for any po > 1, and s > — 1, we have 

l im f£ = M =► lim -4t V ^ P - — ■ (2-2) 

p^oo p s ra— >oo n s+ Z — ' S + 1 



5. If {Xn}^!, {t/ n }5S=i are sequences of real numbers, then for any p > 1 and s,i > —1, we /lave i/iai 

l imsU pM< M , and limsup^<M'^ lim sup V \x p y q \ < M ^ (2.3) 

P-oo P S p^oc p* „-oo 71 po <7^<„ (« + !)(*+!) 

^. {aJnl^Li; {j/n}i^Li ar e sequences of real numbers, andr n is a bounded sequence of positive integer numbers, 
then for all p , qo and s,t > 0, we have that 

Xjy , y v „,. ,. 1 v-^ MM' 
lim = M, and lim ^ = M' => lim — — — > x v y„ = -. 2.4) 

P^oops p^oo pt „^oo nS +t+2 P y « (,s+l)(t+l) V ^ 

Po<P<«, <?o<<2<™ 
|p-9|>r- n 

Before we state the first result of this paper we need to introduce some notations. A path is a string 
A = Ui,.h, ■ ■ -,ji)- A step of A is a pair (j u ,j u+1 ). This is called up if j u+1 > j u + 1, down if j u+1 < j u - 1 and 
a flat if j u+1 = j u . For k > 1, set 

"Pfc = {A = (ji, j 2 , • • • , jfc+i) € Z fc+1 : ji = ife+i, |ju - < 1, 1 < u < k}, 

for the set of paths starting and ending at the same level, and denote by V = Uk>iVk the set of all paths starting 
and ending at the same level. We call A simply a path and we can realize this as a piecewise path taking the 
value j u at u. Now for a given integer p £ Z, we define its shift by p units A + p = (ji + p, j% + P, ■ ■ ■ , jk+i + p) 
and if 1Z is a set of paths in V, we denote lZ + p={\+p:\£ 1Z}. Given a subset fl of Z 2 and a set of numbers 
{«i,j}(i,j)en w e extend this to {a hj } hje z by setting a;j = if (i, j) £ Z 2 \ £1 and a A = a ju j 2 a j2 .j 3 . . .a jk j k+1 . 
Finally, for a given path A £ V , we set max(A) = max{j u , 1 < u < k}, then lj(A) to be the number of crosses 
of the path A with the line y = i + 1/2 and f, (A) the number of flat steps at level i, that is the number of pairs 
ju, ju+i appearing in A with j u = j u+1 = i. For example, A = (-2, -2, -3, -2, -2,-1,0, 1, 1, 0, -1, -2, -1, -2) 
has L_ 3 (A) = 2, L 2 (A) = 4, L_i(A) = 2, 1 (A) = 2, ff_ 2 (A) = 2, ffi(A) = 1, and the other values of If (A), ffj(A) are 
0. Obviously l l+ p(A + p) = and similarly f i+p (A + p) = fj(A). 
Next, define 





= {7 = 


{h,h, ■ • 


■ ■ ,jk,jk+i) £ Vk ■■ max(7) = 0, \j u - j u +i\ = 


1 for 1 < u < k} 




= {7 = 


{h,h, ■ • 


■ ,jk,jk+i) £ V k \ T fe : max(7) = 0} 




r° 
1 fe 


= {7 = 


{h,h, ■ ■ 


■ ,jk,jk+i) £ Vk ■ ji = jk+i = 0, \j u - j u +i\ 


= 1 for 1 < u < k} 


Afc,„ 


= {A = 


Ul,j2, ■ 


■ ■ , jk+i) £V k - I < ju < n, : 1 < u < k} 




A„ 


= U fe >i 


Afe,„ 








= {A = 


{h,k, ■ 


■ -,jk,jk+i) £ ^k,n ■ max(A) = p, \j u - j u+1 \ -- 


= 1, for 1 < u < k} 




= {A = 


Ui,h, ■ 


■ -,jk,jk+i) € A fci „ \ A p kn : max(A) = p}. 





Let's point out a couple of simple properties of these sets. All these paths are strings which move at any given 
step from the previous one by at most one unit and end with the value they started with. Tk is the set of all 
paths of length k with only up or down steps, starting and ending at the same level and staying below the x-axis, 
touching it in at least one point. Similarly, is the set of paths of length k with at least one fiat step, starting 
and ending at the same level and staying below the x-axis all the time but touching it in at least one point. A^^ 
is the sets of all paths of length k staying above the x-axis but below the line y — n. The sets A^ n and A^ n over 
1 < p < n form a partition of the set Ak, n - Notice here an important property which will be exploited below, 
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namely that A — max A e Tfc U T k for any AeP^. In particular for p > k/2 + 1 and any A € we have that 
A — p 6 IV Similarly, for p > k/2 + 1, and any A <E A?'~, one has that A-p€ IT. Therefore, ifp>fc/2 + l, 

AS, n = r*+p and A*- = IT+p. (2.5) 

Consequently, A^ and A^'~ are independent of n for p > fc/2 + 1. This simple property turns out to be an 
important point in proving the next theorem. At last, A„ is the collection of all those paths in V between the 
lines y = and y = n. With these notations, if A n = {cii,j}i<i,j<n is the matrix given in (|1.3|l . then, for any 
path A £ A n , we have 

OA = I II W I I TT % W > (2-6) 




where we use the convention that 0° = 1 . 

In what follows, for a matrix X — (xij)i l j—±,,, n , we denote Tr n (X) = J27=i x a an< ^ = n^- Tn - 
The first result concerns the convergence of the eigenvalue distribution seen at the moment level. 

Theorem 1. Let a > 0. Assume that all random variables d n and b n are independent and there exists a sequence 
{wifc}fc>o, with uiq — 1 so that 



lim E 



(b n /n a f = TO fe for any k>0 (2.7) 



and 

supE 1 ' '" 

n>l 



dn 



< oo for any k > 0. (2.8) 
Denoting X n — -^A n , we have that 

lim E [tr„(X*)] = L k for any k > 0, (2.9) 



and almost surely, 
Moreover, Lk is given by 



lim tr n (X%) = L k for any k > 0. (2.10) 



JO ifk is odd 

Lk = i i rr ti ■ ^- ll > 

{ SS+T ^7er fc IL<o m hh) l f k ls even - 



Proof. Notice that, for k and n fixed, the sets A^ n and A^'~, 1 < p < n are disjoint and Up =1 A^ n UA^'~ = Ak, n - 
As pointed above after the definitions of various T and A sets, for n > p > k/2 + 1, A^'~ = rT + p and 

^fc n ~ Tfe + p which implies that A?'~ and A^ are independent of n. 

Now we denote the elements of the matrix A n by {a,i,j}i<ij< n and then write 



Tr n (A k n ) = Yl 



and since a.ij — for \i — j\ > 1, it follows that 

Tr «(^) = -L E «A, (2-12) 



AGA fc ,„ 



and then 



tr ™( x «) = ^FFrE( E aA + E a * 

p=i V AeA£ n AeA*;" 
We apply Lemma[TJto compute lim^oo E[tr n (X*)]. To this end let's set 

S n= E &A ' S n~ = E &A ' i^-P^ 71 

5 p = E ^ P "= E aA ' ^/2 + l<p. 
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Since k is fixed and 8% — S p ' , S p — S p for k/2 + 1 < p, combined with the fact that ignoring a finite number 
of terms does not change the limit of tr n (X^), we get 



hm E[tr n (A*)] 



1 " \ 

7k+T E (E[^]+E[5^]) 

p=[*/2]+l / 



= 0, 



(2.13) 



where here [sj, stands for the integer part of s. 

Here is the key point of the proof. Invoking (|2.2p and (|2.13|) we reduce the computation of linin^oo E[tr„(J^)] 
to the computation of 

lim —E[S p -} and lim — rELSn 



p— >oo p u 



p" 



To do this, first notice that the sums involved in S p and S p ' are finite, therefore everything reduces to compu- 
tations of the form 

p— >oo 



where 7 is in T k or T k . For the case 7 e r t , according to ((23|) . a 7+p = (]li<o ^j+p^) (lL<o & ! 
products being finite ones. Thus, using the independence of the entries we get 



g\ tl., 



i+p 



\i<0 



Since 7 e r t: , at least one £, (7) is > 1. On the other hand Ylj>o Aj'Ct) + Si<o M7) = ^, from which one gets 
that k - J2ko Hi (7) > 1 an d 



gfc/2 



E[d T+p ] 



< 



1 nf 

p °(fc-£ i<0 li(7)) 



jli (7) 

"j+p 



n E [|(6 s ; +p /^) 1 ' 



(7) 



\i<0 



which together with (|2.7|) and (|2.8p , yields that for 7 e T k 



p— >oo p aK 



(2.14) 



Moreover, since for k odd, Tt = 0, this also shows that 



lim E[tr n (X£)] = 0, 



which is the first part of (|2 . llj) . If k is even, p > k/2 and 7 6 T^, then a 7+p = Hi<o V-H ' an< ^ 



-LEK +p ]=n E [(Wp Q ) M7) 



From this, (|2~7| and ([2~2l . one gets 



i<0 



^ i<0 



(2.15) 



which completes the proof of (|2.1ip . 

For the almost surely convergence, we use Corollary 1.4.9 from |12) . which we state here for reader's conve- 
nience. 

If {X n } n >i is a sequence of independent square integrable random variables and {ro„}™ =1 is a sequence of 
real numbers which increases to 00 as n — > 00, then, for any p > 1, 



V Van 9 pj < 00 =► — V (X p - E[XJ) almost surely. 



p— 1 y 



(2.16) 



P=P0 



From this, it's very easy to deduce the following. 
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Let {X n } n >i be a sequence of square integrable random variables such that there is an integer constant 
q > so that for each r £ {0, 1, . . . q — 1}, {X r+nq \ n >\ is a family of independent random variables. Assume 
also that {w n } n >\ is a sequence of real numbers which increases to oo when n —> oo and has the property that 
lim„^oo Wn+1 = 1. Under these conditions, (|2.16p still holds. 

In our case, we first point out, that almost surely 



lim I tr. 

n — >oo 



1 " \ 

^n)-^5JHT E (S p + S^) =0. (2.17) 

p=[fc/2] + l / 



Let's notice that for each r £ {0, 1, . . . , k — 1}, {S r+nk } n >i and {5 ,r + n ' c . } n >i, are sequences of independent 
random variables. Now we take w n = n +1 . We show first that 

^ var(S*P) , ^ var (SP-) 

1^ -pl^kW < °°' aild 1^ p 2 a k+2 < °°> 
p—k p=k 



which follows once we know that for any 7 £ Tk U T 



k • 



p—k 

To prove this, from ([21]), a 2 +p = {Uj< af^) (llt<o tf+p^) and usin S <E3) and (Ml ° ne S ets that for a 
certain constant Ck > 0, 

^E[a* +p ] <C- fc forp> fc/2, 



P 

which is enough to justify (*) and thus, by (|2.16[) . that 



1 - 

^ £ (S>-E[i?\+S»--E[S>--])-^0. 

p=[fc/2]+l 

This together with (f2~T7|) . ([2~T4|) and (f2TT5|) prove (|2~T0ll . □ 



Remark 1. Condition (|2.8p can 6e relaxed under the assumption that are the moments of a measure which 
satisfy Carlemans condition X^fcLo ^k ^ < 00. In this case, if we replace the condition (|2.8[) by the condition 
that 

supE[|d„|] < 00, 

n>l 

we can conclude that the empirical distribution of eigenvalues of X n converges to the measure whose moments 
are given by . 

The proof of this fact is basically given in J^j, page 615 where it is proved that the diagonal part can be 
removed. The only essential fact which is needed there is that sirp p>1 P(\d p \ > n a e) — o(l) which follows from 
the above condition and Chebyshev's inequality. 

Corollary 1. Within the notations of the theorem above, assume that rrik = 1 for any even k. Then the numbers 
Lfc are the moments of Ullman's distribution v a {dx) = h a (x)dx with (a > 0) 

1 f 2 t~ 1+1 f a 
h a (x) = 1[_2,2](V) — / -j==dt. 

a7r J\x\/2 V4 - t 2 

These are obtained as distribution of ' T a W , where W has the arsine law distribution (l[-2,2] ( x ) i2 rfjJ and 
T is an independent uniform on [0, 1]. 

In some cases, closed formulae are available, as for example, 

[M-2,2]{x) {2+x2 lf^ a = 1/4 

h a (x) = I ±.l [ _ 2a] {x)y/j—& a = 1/2 (2.18) 

[^l[- 2 , 2] log((2 + V4^)/H) a = l, 

In particular for the model (|1.2|) . the limiting distribution is the semicircular law. 
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Proof. Since rrik = 1 for all k even, we have that the products in (|2.1ip involving m's equal one. The number of 
such terms is given by the number of paths in r&, which turns out to be (Jj 2 ) ^ or ^ even - One very quick way 
to see this is that any path 7 in is perfectly determined by the prescription of the places where the up steps 
start, the rest of the positions being filled in with down steps. Since the path must have the same starting and 
ending point, it means that there are exactly k/2 up steps. The way of choosing k/2 positions out of k points is 
just ( fc y 2 ) • This means that 

f k odd 

Lk = <( M (*) 
1%T Severn 



For k even and a — 1/2 these are the moments of the celebrated semicircle law ^:l[-2.2]( a; )\/4 — x 2 dx. 

Even though the semicircle plays and important role here, it is the case a = which is the most important 
one. For a — 0, we have Lk — ( fe / 2 )- One can cnec k directly that the measure having these properties is 
the measure i>o(dx) = 1[— 2,2] ( x ) x ' } < ^ x ' Now if ^ is a random variable with distribution Vq and T is an 
independent and uniform on [0, 1], then T a W has the moments given by (*). From here the rest follows by 
direct calculations. □ 

Remark 2. The system (|2. 1 1|) is invertible in the sense that for any given the sequence Lk, k even, one can 
solve uniquely for the sequence rrik, k even, since the system (|2.11| is a triangular one. To simplify the notations, 
set Mk = (ak + l)£fc. Then we have for the first lines of the system (|2.11[) : 

M = m = 1 

M 2 = 2m 2 

M A = 2m 4 + Am\ 

Me = 2ma + Ylm^m\ + §m\ 

Ms = 2ms + 1677167712 + 127774 + 3277147772 + 8t772- 

We can solve for m 's in terms of M 's in this case as: 

m = M Q = 1 

777 2 = \ M 1 

m 4 = i(M 4 -M 2 2 ) 

77i 6 = - (4M 6 - I2M4M1 + 9M|) 
8 

77i 8 = i(2M 8 - 8M 6 M 2 - 6M 4 2 + 28M 4 M 2 2 - 17M 2 4 ). 

It is of interest a combinatorial interpretation of this inversion. Moreover, one such interpretation could perhaps 
lead to an analytic interpretation, one which would allow characterization of the situation in which the numbers 
Lk are the moments of a real measure. 
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Next, we would like to compute the limiting distribution in one particular case in which b n /n a converges, 
not to a constant, but to a Bernoulli random variable. The next proposition also shows that the numbers Lfc are 
true moments of a measure under some reasonable conditions. 

Before we state the main result we introduce a class of infinite random matrices known somehow in the 
theory of random operators as the Anderson model ([13j). Assume that {X n } ne z are given bounded iid random 
variables. Then we define 



A = 




(2.19) 



where the marked element is the (0, 0) element. We can realize this matrix as a symmetric random Jacobi 
operator acting on £ 2 (Z). 

Consider the unitary map from U : £ 2 ("E) — ► L 2 {S V ), where S 1 = {z E C : \z\ = 1} is endowed with the 
uniform measure and U({di}i£z)(z) — 'Yl,i^j J a i z Then the matrix A becomes the random operator which is 
given by A = UAU' 1 

A{z l ) = {XiZ' 1 + X t+1 z)z\ (2.20) 
Finally if e, € £ 2 (Z) is the vector with 1 on the ith component and otherwise, then 

Ai.o = M- e 0) e o)^ 2 (z) = {A l,l) L 2 (s i), 
where here 1 is the constant function 1 on S 1 . 

Proposition 1. For a > 0, assume that there is a bounded random variable Y with moments K[Y 2k ] = m^u for 
k > 0. Then there is a bounded random variable Z whose kth moments are Lk- 

The distribution of Z is of the form v{dx) = rf(x)dx + (1 — t)5q, where < r < 1 and f is a measurable 
density function. 

Ifa>Q,0<9<l and mik = & for all k > 1, then the distribution whose moments are Lk, k > is given 

by 

'8 (dx) 9 = 

v e , a {dx) = { v a (dx) = 1 (2.21) 



where 



and for M > 3 



fe,a( x ) 



H-2,2) 



M- g f e ^(x)dx + ^§6 {dx) < 9 < 1, 
(l+0)(l_0) a ^ 



(x)- 



N>1 



g N+2 {x), 



p/arccos(|a;|/2)/7r] 
g M ( X ) = la-jVa-l £ 



^2 1 / a \cos(uTr/M)\ 1 / c 



with the convention that = and [s] stands for the largest integer < s. 



Proof. Assume that the distribution of Y is a measure fi with support in the closed finite interval /. Consider 
now the probability space fl = I z and P = the product probability on O. We denote by uii the ith 

component of to. Then, define the Hilbert space H = {x = {xj}i S z : Xj € L 2 (Q,, P), X^ez ll 3 ^' Hl 2 (£"2 p) < °°} 
with the scalar product given by (x, y) = J2jez ^ P [ x jVj]- ^n this Hilbert space we consider the operator A 
given by 

Ax = {MjXj-i + M j+ iXj + i} je z 

for x = {xj}j£%. Here : L 2 (fi, P) — > L 2 (£l, P) is the multiplication operator given by (M.ix)(uj) = ujix(ui) for 
any x S L 2 (fl, P). Since / is a closed finite interval, M is a bounded operator and this in turn yields that the 
operator A is also a bounded selfadjoint operator. 
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1 2-2-1 



2 -2 -1 



( a ) /o.4,l/4 (b) /o.6,l/4 ( c ) fo. 95, 1/4 

Figure 2: Density /e,i/4 for # = 0.4,0.6,0.95 and the smooth density h\u. 






(a) /o.4,l/2 



(b) fo. 6,1/2 



(c) /o. 



95,1/2 



Figure 3: Density for = 0.4,0.6,0.95 and the smooth semicircular density h\ii = ^:\/4~ 




(a) /0.4,1 (b) /o.6,l (c) /o.95,l 

Figure 4: Density f@i, 9 = 0.4,0.6,0.95 with the smooth density hi 



Now we define e = {ej}j£%, where ej = 1 if j = and ej — otherwise. We will prove that 

'0 k odd 



<A fc e,e) = 



,E 7e r fc n i <o m i I (7) keven 



(2.22) 



To do this we first take the random variables Xi : Cl — ► K given by -Xj(w) = Wj. The set {JQ}j 6 z is a set of iid 
random variables with distribution \x. With the random infinite matrix A given by (|2.19[) . notice that 



(A*e,e)=E[^„] 



(2.23) 



which means that we first compute formally A k and then take expectation of the (0, 0) component. From this, 
if we use <Zjj for the entry of A, then 



A 



0.0 



■ a i k .O 



fcodd 
2~] 7 er° a i k even - 



This, together with the fact that a^j+i = -Xj+i, = that {JQ}j 6 z are iid with distribution and a 

moment of thinking, gives for k even, 



E Ko] = E I1 E [^ (7) ] = E EK(7>> 

7 er°iez 7 er fc i<o 



(2.24) 
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which proves (|2.22[) . 

On the other hand, since A is a bounded selfadjoint operator, we can take its spectral measure £,(dt) and 
then w(dt) = (£(dt)e,e). We then have that J R t k xu(dt) = (A k e,e). Now if we take a random variable W with 
distribution vj and T an independent uniform random variable on [0, 1], one can check, using (|2.22|) and (|2.1ip . 
that Z = T a W has the moments Lfc. It is an easy exercise to verify that the general form of such distributions 
is v(dx) = rf(x)dx + (1 — t)Sq with r £ [0, 1]. 

For the second part, the case 9 = is obvious. Even though the case 9 = 1 is covered by Corollary [TJ we 
want to employ the arguments used in this proof to reprove it. The random variable Y in this case is simply 
the constant 1. Therefore, the operator A becomes a nonrandom operator and, using the representation given 
by (|2.20p . is in fact the multiplication with z + z^ 1 on L 2 (S' 1 ). Consequently the spectrum is [—2, 2] and 



x k p(dx) = (A l,l) L 2 (s i ) 



-i\fc 



dz 



1 

27 



2 fc (cos s) k ds 



1 



- T dx, 



IS 1 *' n JO 71 J -2 \f£— X- 

which results with w{dx) = vq = ^l[_a,2] ( x ) v / 4 1 p given in Corrolary [1] . From this a simple calculation shows 
that the distribution of T a W where W has distribution w is the v a given in Corrolary [TJ 

Next, if 7712/t = 9 for all k > 1, it is easy to see that Y whose even moments are uiik is a Bernoulli random 
variable with probability 9 of 1 and probability 1 — 9 of 0. Thus, the matrix A has elements or 1. To compute 
the (0,0) entry of A k , we notice first that (observing when the first appears in the sequences {Xi}i<o and 
{Xi}i> ) 

l,m>Q 

where *4;, m is the matrix A, with Xq = X_i = ■ ■ ■ = X- m = 1, X_ m _ 
Xi + i = 0. Thus, the matrix Ai, m is a block matrix of the form 



and Xx = X 2 



X, 



1, 



A. 



l.m 



c, 










Br 






D 



with the square matrix 



1 



1 1 
..10 



the marked entry being the (0, 0) entry of the matrix Ai >m and counting from it, there are I + 1 rows to the 
top and m rows to the bottom. Bi, m is a m + I + 1 tridiagonal matrix with only 1 on the upper and lower 
sub-diagonals and otherwise. The key point is the fact that for N > 0, 



f^ N B k Nl 



l,m>0 
l+m=N 



where Bn is any of the (TV + 1) x (N + 1) matrices Bi lTn with I 
other hand, 

N 



i = 9 N (l-9f E (^0,0 = ^(1 

m — N and the marked entry removed. On the 



i,m>0 
l+m=N 



Tr N B h N = E 



n k 



u=l 



where 7] Ui n are the eigenvalues counted with their multiplicity. Now, if qN(x) — det^J/v — Bn) is the character- 
istic polynomial of Bn, an easy induction argument shows that 

qn(x) = xq N -i{x) - qn-2{x), 

with qo(x) — x and q\(x) = x 2 — 1. These, up to scaling, are the Chebyshev polynomials of the second kind. 
Precisely, we have qN{x) — Un+\{% )■ As it's well known, the roots of Un are cos (^jf^fij, 1 < u < N and this 
shows that the eigenvalues of Bn are 

rj u ,N = 2 cos 



N + 2 



fori < u < N + 1. 
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From here one gets that 

JV+l 



JV>0 u=l 

where 

N+l 
AT>0 u=l 

If this is the distribution of W and T is uniform on [0, 1), then T a W has the distribution given by 
1—6 N+1 \x\ 1 / a ~ 1 9 N {1 — 9) 2 

U e>a (dx) = ——S (dx) E nOl/ct+llmf ™ ^1 / cos(u7r / (N+2))\,2\ cos(u7r / (N+2))\)(x)dx 

0+1 N>1 u=l \ coi >{W+2)\ 

and from here, rearrangements bring this to the form given in (|2.21|) . □ 

Theorem [1] gives the zero order convergence in moments of the distribution of eigenvalues. Here we are 
interested in the first order convergence. The statement can be made more general, but for the sake of simplicity, 
we give the next theorem in the following form. 

Theorem 2. Suppose that the conditions (|2.7|) and (|2.8|) from Theorem]^ hold. In addition, assume that there 
exist < v < min(l, 2a), and numbers ad and £fc so that 

lim n v (E[(b n /n a ) k ] - m k ) = £ fe , k > 0, (2.25) 

n—*oo 

E[d n ] = 0, Vn, and lim var[d„] = a 2 d (2.26) 

n — >oo 

then, linin^oo n v (E[tr n (X^)] — Lk) is given by 

' ak -v+i E 7S r fc E i<0 6/(7) IIi<o,iyy m )U(7)> fc eue " > 2 ^< mm (!' 2 «) 

S^+lWE 7 gr»'- Ili<o m ii(7) +E 7 er fc E J <o?i J (7)n i <o, I ^ TO i,(7))' fc even >2,v = 2a<l 

^T^Lk + I E 7 er fc (E*<o *1*(7)) ELko "^(7) + 3 E 7 er fc E J<0 6,(7) IL<cm^ m Ui)> k even >2,v = K2a 
^1T-L k + ±E 7 er fc (Ei<oJ' 1 i(7)) IL<o m ii(T) 

E 7er ='- IL<o TO i,(7) + f E 7 er fc E,-<o 6,(7) Ili<o,i^ m M7)' k even >2,v = l = 2a 

0, k odd, or k = 0, 



where IV is i/ie se£ 0/ paths in T k with exactly two flat steps, both at the same level. 
Proof. We use the notations from Theorem [TJ Notice first that, 



(2.27) 



\ p=[fc/2] + l / 



= 0. 



Thus, to prove (|2.27|) it suffices to find 

\ \p=[fc/2] + l 

and, according to Lemma [T] and the definition of Lfe, it reduces to 



lim -^J E [a 7 +p] for 7 G T fc , (*) 



and 



/„afc+l _ f„ _ 1 \afc+l\ \ 

SS,=5S=5- *K +P ]-^ ^ l U m H,) > fOT Ter, (**) 



i<0 
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Finally this can be done by using 
1 



If 7 € T. , then, 



V 



J<0 



.(7) 



\i<0 



If 7 has exactly one flat step, this quantity is zero because E[d n ] = for all n. On the other hand if 7 has exactly 
two fiat steps, then for v < 2a this whole term goes to with p. In the case 7 has more than 3 flat steps, one 
gets that J2j<o %(7) — 3 and so again the term goes to 0. The only case we get something nonzero is the case 
when v — 2a and 7 has exactly two flat steps at the same level. In this case, according to (|2.26p . one has that 



p— *oo 



i<0 



Summing over all possible paths in T k 

1 



lim 

p — >oo p° 



[s p -] = 



°"dE 7 er 2 fe - ]li<o m ii(7)' v = 2a 



otherwise. 



If 7 € r fc , then, 
' ' Ef. 



t 7+PJ 



ctfc+1 



(p-1 



afc + 1 



n 

i<0 



m i;(7; 



^ V<0 



(7) 



A' 



P "II mi >(7) 

i<0 / 



( p ak+l _ ( p _ !)afe+l _ + l) p <*kj 



II TO 1.(7) 



(ak + l)p ak ~ v 

^ ( n E [(w(p+o q ) i,i(7) ] ) (p v (e [(6 l+P /P a ) iiW 

i<0 \/i<i / 

+ E f II E [( & ^+p/(p + *n lh(7) ] I (p w (a + */p) Ql,(7) - 1 

i<0 \h<i / 



m li(~f) 



ri TO i S (7) 




(p at+1 - (p - 1) Q ' £ + 1 - (afc + l)p Qfe ) 
(afe + l)p Qfe - u 



i<0 



> [£j<0& 3 -(7)IL<0,i#j m l i (7)> W < 1 

P -oo I +aY /i<0 *U(i)) Uj<o m hh) + Ej<0 6,( 7 ) II, „...., "»1 4 ( 7 )> v = L 

The rest of (|2~2"Tj) follows. □ 

Remark 3. We mention here the following equalities for k even k > 2, 

|rM = fc2 fe - 3 , 

E l ^) 2 = fc2 fc , (*) 
7er fc 

£ ^il i ( 7 ) = -fc2 fe - 1 . 
7er fc i<0 

Applying these to the model from (|1.2|) . we #e£ £/ia£ lim^oo n(E[tr n (X^)] - fe/2+1 (fc/2) ) = (f — f)(2 fe_1 — ( fc / 2 )) 
which are the moments of the measure \{5-2 + ^2) — %— 2,2] ( x ) i2 rfi- This is J21 Lemma 2.20] and we just 
reproved it. 

The proof of the first equation can be done by counting all the paths by directly. The second equation can be 
proved using the model Alt) ([2~T9]) with X, = 1 + tYi wh ere Y% is a sequence of iid N(0, 1) and use (|2.24j) to get 
that the identity we are looking at is just the coefficient oft 2 of E[(A(t)) 0() ]. Then we can write A{t) = B + tC. 
where B is the matrix with 1 on the subdiagonals and C has iid normal N(0, 1) on the subdiagonals. One can 
compute the coefficient of t 2 in the kth power of A(t) as a product of the form B kl CB k2 CBk 3 . The powers B 



15 



can be explicitly computed and then the rest is simple combinatorics. The proof of the third equality in (*) can 
be done in the following way. First realize that the term ^2 i<0 Hi (7) is the negative area between the path and 
the x-axis. Then one can decompose any path in Tk as two Dyke paths with certain properties. Finally, one can 
count the number of paths with a certain area (see J3j Proposition 6]), together with manipulations of generating 
functions to get the equality in there. 

The outlined proofs are long and ad-hoc. More direct and natural combinatorial proofs are desirable though. 



3 Fluctuations 

Under the conditions in Theorem [1] we have almost surely the convergence of the distribution of the eigenvalues 
of X n . In this section we are interested in the "fluctuations" from the limiting distribution. Theorem [TJ states 
that almost surely, 

lim (ix n X* - E[tr B X*]) = 0. 

n — >oo 

Next wc arc interested in how this happens. More precisely, what is the right factor we should multiply tr n X„ — 
E[tr„X^] with to make this converge to something? Assume that we multiply this by n v with n > 0. What is 
going to be the right 77? Let's take a look at the case k = 2. Then, 

n n— 1 

n" (tr n X 2 - E[tr n X 2 n }) = £ [d] - E$]) + £ (b 2 - E[b*}) . 

i=i j=l 

Now, for any < n < 2a, the first sum of this goes to by the Strong Law of Large Numbers. The other sum 
can be written as 



^ E ^ - m) = E ^ 

3=1 " 3=1 



with Yj = b 2 /j 2a - E[b 2 /j 2a }. 

Let's assume that {Yj}'^ 1 is a sequence of independent random variables so that in distribution sense j e Yj ~ 
U for some a > e > and U a zero mean random variable with variance var(C7) > 0. Then we are looking at 
the condition that 

-. n— 1 -. n— 1 

7 - 1 \" Jay ^ 1 Ja-err 

Zjn ~ n l+2a- V J '■> n l+2a-r, U 3 

3=1 3=1 

is converging in distribution (here {Uj} are iid with the same distribution as U). Take ip{x) so that the charac- 
teristic function of U is E[e ltu ] = e^*> with ^(0) = and ip'(0) = 0. The condition of convergence is translated 
roughly as convergence when n — > 00 of 

n-l 

^(tj 2a - e /n 1+2a -v). 
j=l 

Since U is not constant 0, this implies that ^"(O) = var([/) > 0. Now, Taylor expansion ip(x) = ip"(0)x 2 /2 + 
0(x 3 ) yields 



n-l 1 n-l ( var(t/)t 2 _ -, /t> 

^^(^-7n 1+2 ^)^var(^_i_^,4- 272r j „ <l - 



i=i 3=1 



2(l+4a-2e) 

or 00 otherwise. 



Therefore the choice in this case is obviously n = 1/2 + e. Moreover, this also shows that the limiting distribution 
of Z n is normal. 

Another way of guessing 77 is from the general statements of CLT, for variables which are not necessarily 
identical. 

Before we state the next result, we need some definitions. 

We say that the paths Ai,A 2 S V do not share a level if for any i e Z, li(Ai) 7^ implies li(A 2 ) = 0, 
MA2) 7^ implies li(Ai) = and similarly for the flat levels, ff;(Ai) 7^ implies ffi(A 2 ) = and ffi(A 2 ) 7^ implies 
fj(Ai) = 0. We say that Ai and A 2 share a level if there is an i so that both li(Ai) are lj(A 2 ) not zero or both 
ffj(Ai) and fi(A 2 ) are not zero. 
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For k, I > 1, we set r(fc, I) by 

{{(71 j 72) G T-'fc x Pi '■ max(max(7i), max(7 2 )) = 0, 71, 72 share a level and have no flat steps}, if k, I even 

{(7i) 72) e x : max(max(7 1 ), max(7 2 )) = 0, 71, 72 have exactly one flat step each and is shared}, if k, I odd. 

(3.1) 

Remark 4. For k,l odd, the number of paths in T(k,l) is the number of paths of (71,72) G Pk x Pi, 7i>72 
having exactly one flat step on the x-axis. These pairs can be constructed as follows. Pick two paths 7^ and 7 2 
of length k — 1 and I — 1 with only up or down steps. Then insert any flat step at any level and move the paths 

so that the level steps are on the x-axis. For j[, there are choices for the path and k ways of inserting the 

2 

flat step. Similarly for 7A, so the total number of paths in T(k,l) is kl( h~i ) ( 1-1 ) ■ 

2 2 

Theorem 3. In addition to conditions (|2.7[) and (|2.8[) of Theorem^ assume that, 

lim var(d„) = a% (3.2) 

n — >oo 

and there exists < e < a so that if z k = n e {ib n /n a ) k — K[b n /n a ) k ] s j , then for k, I > &o£/i even, £/iere exists 
C(k, I) such that 

lim cov(z£,z^) = C(fc,Z). (3.3) 

— S-OO 

and 

su P E[|z£n < oo,Vfc,m> 1. (3.4) 

n>l 



./Vow se< 

' f = I) : = 



' p = n . = J a ( fc +z)+i E( 7l , 72 ) S r(M) (rife<o m i fl (7i)+U(72) - rih.<o TO U(7i) TO U(72)) . Mfc, ' even 
0, otherwise. 

"'2 + ' y v> g(l.(7i),l,(72)) fork! ever, 

a(fe+i)+l-2e Z^( 7l , 72 )er(fc,Z) Z^i<0 lh(7i)+>h(T3) ' ^ eye,t 



D(fc,0 



< e < a 



0, otherwise. 

V V C{h(li), 1,(72)) jr i, 7 

Q (fe+;) Z^( 7ll72 )er(fe,0 Z^i<o m J h (7i)+i h (T2) ' J"' K ' 1 euea 



k + l 

m 



J d" L 2 



fc_i) (i-i), fork, I odd, 



*(k+i) 

otherwise. 

(3.5) 

For any polynomial P(x) = wq + wix + • ■ - + wj^x N denote S n {P) — n e+1 / 2 (ti n (P(X n )) — E[tr n (P(X n ))]) . Then 

lim S„(P) = 7V(0,a(P) 2 ) 

n — »oo 



<t(P) 2 = J2 w k wiD{k,l). 



k,l>l 

In particular, if S n (k) = S n (x ), this implies that the family {S n (k)}k>i converges in moments to a Gaussian 
family {S(k)}k>i with covariance function D(k,l) and 

lim S n {k) = N(0,oi) (3.6) 

n — >oo 

where a\ = D(k, k). 

Remark 5. Let's point out that in the case < e < a, condition (|3.3[) implies in particular that mk+i = mkmi 
and hence in probability too. 



for any k, I even. This in turn means that = m^ 2 for any k even, or that b^/n 2 " ► m 2 in distribution 



Proof. Write 

N 

(S n (P)) j = ^ w kl w k2 ■ ■ ■Wk j S n (k 1 )S n (k 2 ) ■ ■ ■ S n (kj). (3.7) 
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Since there is a finite number of terms in the above sum, to study the behavior of E (S n (P)) 
deal with 

hm E[S n (k 1 )S n (k 2 )...S n (k J )} 



, it suffices to 



for a given sequence k\ , k% 1 ■ 



, kj . Now, since 



S n (k) = 



n (afe+(l-2e)/2) 



(o A -E[a A ]) 



we have 

E[S n (k 1 )S n (k 2 )...S n (k J )} 



l23=i(°*i+(l-2e)/2) 



E [(oa, - E[oA 1 ])(a A2 - E[a A2 ]) . . . (a Aj - E[a Xj ])] 



A;GA fci ,, 
l<i< 3 



(3-8) 

Next we define a notion of connectedness for paths. A set of paths C is called connected if for any two paths A 
and A' in C there are paths A = Ai, A2, . . . , X u = A' in C so that Aj, Aj+i share at least a level. Otherwise we say 
that C is not connected or simply disconnected. The notion of connectedness in this context is an equivalence 
relation. Therefore any set C can be written as a disjoint union C\ U C 2 • ■ ■ U C p where each Cj is connected. The 
sets Ci,C 2 , ■ ■ ■ ,C P are called the connected components of C. If A is a path in C, then the connected component 
containing it is the set of all paths which can be connected with it. In particular any two paths from different 
components do not share a level. 

With this concept at hand, we return to (|3.8|) and split the sum in sums over all connected components. Then 
we organize the connected components in the following way. For a given partition A of {1, 2, ... we consider 
C\ the subset of (Ai, X 2 , . . . with Ai £ ^k i: n so that Cg — {A.; : i S 8} are the connected components of 
where 8 runs over all elements of A. In other words the connected components of (Ai, X 2 , . . . , Xj) are indexed by 
the subsets 8 £ A. Now, since any two connected components are disjoint, combined with the independence of 
the entries of the matrix A ni justifies the following rewriting 



E[S n (k 1 )S n (k 2 )...S n (k j )} = 



1 



i£i =I («*i+(l-2e)/2) 



Next we fix a partition A of {1,2, 



A ( Al , A2 ,..., Aj -)eC2, <5eA 
,j}. The idea is to find the limit of 



]J (o A — E[<2 A ]) 
Aec™ 



U A = 



1 



>a=i(«fe+(i-^)/ 2 ) 



E IP 

( A i, A2 ,...,A 3 )ec?s <5eA 



n («a - e[o A ] 

Aec 6 - 



To clarify and explain the main idea let's introduce first some notations. For a given k = (fci, k 2 , . . . , kj) and A 
a partition of {1, 2, . . . , j}, we set 

r(k : A) = {T = (71, 72, . . . , 7,-) : j u € Vk u , and for 8 € A, max(max(7„)) = 0, : u € 8} is a connected set}. 

This is the set of all paths under the x-axis so that by isolating the paths indexed by 8, we obtain a connected 
set with the maximum of all heights being 0. 

Notice that for a given A, the set T(k : A) is actually in a one-to-one correspondence with the set x^gAlXka : 
<5), where k,5 is the vector k with the components which do not belong to 8 removed. Obviously there is a finite 
number of elements in T(k : A). Now if we take a connected component in Cf, with max ue ,5(max(A u )) = p$, 
then ((X u — ps)ues)seA ^ x <5eAr(k,5 : 8). Ignoring eventually a finite number of terms in the expression of U^, 
the limit of U A is the same as the limit of 



V A = 



1 



rer(k:A 



E IP 



E[, 



]) 



where the set fl n (r) = {(ps)seA ■ n > ps > Po, so that {75 + p s } n {75- + ps>} = 0, for any 8 ^ 8'}, p Q being 
a fixed large number depending only on 7a and k. If the set A has just one element, namely the whole set 
{1, 2, . . . ,j}, then the above sum takes the simpler form 



V. 



A 



E 



E^ 



n E„ e{1>2 3} (afe u + (l- 2e )/2) 

rer(k:A) n p=po 



u£{l,2,...,j} 
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Using Lemma [TJ wc can find this limit once one can compute the following 



W(T) = lim 



p (E u6 {i,2,...,j}(«fc«.+(l-2e)/2))-l 



«e{i,2,...,j} 



(3.9) 



Once we know this, we can go back to the case of an arbitrary partition A and use part 3 and 4 of Lemma [T] to 
compute the limit of V„. Here are the formulae. For A with just one component {1, 2, . . . ,j}, we have 



lim V A 



E 



2W(T) 



rer(k:A 



2a|k|+j(l-2e) 



with |k| = fci + ki + ■ ■ ■ + kj. The general formula which follows from this and a repeated application of Lemma[T] 
is that for an arbitrary partition A we have 



i^„ a = e n 



2W( 74 ) 



rer(k:A) <5eA 



2a|k 5 | + |<y|(l-2e) 



(3.10) 



Now we want to compute W^(r) when T e F(k : A) and A = {{1, 2, . . . ,j}}. In the following, for a set we 
denote the number of its elements by |f2|. 

Case 1: j = 1. In this case due to the fact that E[(a^ — E[a>])] = 0, we get that U£ — and in 
particular, W(T) = also. 

Case 2: j > 2. We show in this case that W(T) = 0. 

To do this we will prove something more general. Namely we show that for a fixed F 6 T(k : A), 



E 



II ( a 7u+p -E[a 7u+P ]) 
«e{i,2,...,j} 



Oip 



(3.11) 



For any path 7, recall (|2.6|) which takes the form a 7+p = flI fl <o ^3+?) (ll/i<o ' tn ^ s P rocm ct being 

actually a finite one. To make the writing in a reasonable form for the expansion of the left hand side in (|3.11|) . 
we rewrite 



(7) 



where 



oil) 



d i < 



b i>0 

Since the entries are independent, we have 



a 1+p - E[a 7+p ] = II d 



■U~f) (Ah) 
g+p \ a g+i 



and Txii 



— E 



fi i<0 
L, i > 0. 




f 9 (7) 



i<0 \a<0 / \/i<i 




tlh(7) 



which can be rewritten as 



%7+p] =eii ( c a m -i?) T,(i) ( c (cf? - e he!?]) 1 "" ( e K m -N } ]) 



X 7+P 



zez iez 



, (*) 



with 



Jl Kl \l l = l \l i>l 

n{l) = i>/' ^ W= ow= i<i. 



Notice that 77 + vi + £j = 1, X^ez^W = 1 f° r am/ ' e ^ an< ^ if 7 £ T^k, then Xiez mi (7) = ^- Using these 
formulae for 71,72, ■■■ ,7j, after multiplying out the factors, the left hand side in (|3.11[) becomes 



e 

ii,J2,...,iiGZiGZ 



n 



nii(7„) 
■H\ 
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Now, if i < 0, then c(i) = d and then (|2.8p combined with Holder's inequality yields that each product in the 
above sum with i < is bounded by a constant. If i > 0, then c(i) = 6, in^ = l_i and in this case, using Holder's 
inequality, (|2.T[> . (|3.4[) and (|3.3|) one can show that 



E 



i(7u) 



< (7pEi=i( am i(7 u )-ei / i„(i)) j 



where the constant C depends only on the paths 71, 72, . . . , 7/. This means that for fixed li, I2, ■ ■ ■ , lj € Z 



n 0) 



mij(7«) 
P— 1*1 



T l u (i) 



nii(7„) 



E 






(e 




) 




c Wd-|«| 






C Wp-|i| 





,(l„(i) 



< Cp52l=i E,>o( Qm .(7»)-fi„('))^ 



Next we have 



afc„ — e if 7„ contains no flat step 
(airii(7 n ) — ei/; u (i)) < < afc„ — a if 7„ contains exactly one flat step 

ak u — 2a if 7„ contains two or more flat steps. 



i>0 



ue{i,2,...,j} 



< { 



To see this, one should notice that if 7„ does not contain a flat step, then X^>o m i(7ii) = 2~^i<oM7'«) = 
while ^2 i>0 ^i u {i) = 1- In the case 7„ has just one flat step then, if l u < 0, then ^2 i>0 h(ju) = k u — 1 while 
Si>o (*) = ^ an( i ^ > 0) then 2~Zi>o ^(Tu) — ^« — 1 while X)i>o ^ W = 1 which justifies the first part. In 
the case 7„ has more than one flat step, then J2i>o l»(7u) < ^« — 2 and the rest follows. Hence, since e < a, 

CpI?L=i( ak »- e ) if any7„, 1 < u < j, has no flat steps 

Cp(Ei=i(«fc«.-s))-(«-e) if any7M; 1 < u < j, has at most one flat step 

and at least one has exactly one flat step 
(jp\jy^i{ a k u -e))-a •£ Qne f 1 <u< j, has two or more flat steps. 

(3.12) 

which suffices to prove (|3.11[) . 

Case 3: j = 2. In this case, T = (71,72) and we need to compute 

p 1 ™, pa fc 1+ ifc 2 -2 £ E [( Q 7i+p - E Ki+p])( a 72+ P - E[a 73+P ])]. (3.13) 

Here we distinguish the cases e < a and e — a. From equation (|3.12[) . we see that for e < a, the dominant 
term is the one involving only sums over the paths with no flat steps. If e = a, then we need to consider also 
the paths with exactly one flat step. 

First we consider the contribution from the paths with no flat steps. To carry this out, invoke (|2.6p and since 
there are no flat steps, a 7+p = Yl h<0 ^h+p an< ^ 

a 7+p - E[a 7+P ] = E (il ) ifiS 

from which one gets 

( a 7l+p~^[ a 7l+p])( a 72+P — E [ a 



i<0 \h<i 




e nc i)+ih(72) 

n<»2<0 \h<ii 



x n 

\ii <h<i 2 



'lh(7i 
^+P 



i<0 \7i<i 



lh(7l)+lft(72) 
>i+p 



klh(72) 

°7l+p 



f 6 l*(7l) 



-E 



2(72 



M11) 

"i+P 



E 



.1.2(72) 

»2+P 



E 



lh(7) 
h+p 



,ln(7i) 
ii+P 



JM72A 
n+p J 



(1(72) 
/i+P 



E 



J.lh(7i) 
°h+p 



)(* 



i 72) - E 
+p 



\i2<h 



M12) 
U i+P 



n4» 



(72) 
h+p 



E 



■lh(7i 
'h+p 



\i<h 



e nc +uhi) ((« 



72) 



E 



»2<il<0 \h<i2 

>< ( n e 



.1. 2 (72) 
'12+P 



MM 

h+P 



'lfc(72 



1.1/. (7i 
^+p 



^1(71) 



E 



.lii(7i) 
°ii+P 



f.(7l) 
h+p 



\ii<h 



E 



■lfc(72) 

; /l+p 



20 



After taking expectation in this formula, from the independence of the entries, one arrives at 



E 



ilh(7i) 
°h+p 



cov( a71+p , a72+p ) = £ n E [C l)+Ift(72) ] -v(^\^) J[E [ b {f; 

i<0 \h<i I \i<h / 

from which, according to (|3.3p , it follows that 

j™, Da fc 1+ ifc 2 -2 e C0V ( a 7i+P' a 72+p) = I II m U(7l)+l fl (72) J C(li(7l) 5 li(72)) I n m ^(7i) m U(72) ) ■ ( 3 ' 14 ) 

~* ^ i<0 \/i<i / \i<h / 

Let's point out that in the case e = 0, one has C(fc, I) = rrik+i — rrikmi and in this case the formula simplifies to 
the one given in (|3.5p . 

Next we deal with the case in which there are fiat steps in 71 and/or 72. Thus we need only consider the 
case e = a. In the first place if only one of them has a flat step then we may assume that 71 has one flat 
step and 72 does not. Then we write a 7l + p = d g+p Ylh<a ^h+p where g is the level of the flat step. Hence 



a 



Ylh<o an d then because of the independence of the entries, 



72+P — li/i<0 u h+p 



cov (a 7l + p , a 72 + p ) — E[d s +p]cov(a 7 j +p , a 72 + p ), 



where "f[ is the path obtained from 71 by removing the flat step and gluing together the remaining parts. As a 
consequence of the above, the path (71,72) does not have a flat step and for the T = (71,72) we can use (|3.14p . 
Taking the limit over p — > 00 and keeping in mind that now the path j[ has length k\ — 1, using the previous 
step we get that 

p 1 ™ pa (k 1+ k 2 -2) COV Ki+P' fl 72+ P ) = 

The third situation here is the one in which both 71 and 72 have a flat step. In this case we write a 7l+p = 

a gi+p a 7l+P — d 92 +P Hh<o & l+p 2) where 31 and g 2 are the levels of the flat steps in 71 and 72. 

Denote by j[ and j' 2 the paths obtained by removing the flat steps. In this case T — (7i,7 2 ) has no flat steps 
and the length of j[ is k\ — 1, while the length of 7 2 is hi — 1. Now a simple calculation gives 



cov 



(a 7l+p , a 72+p ) — E[d gi+ p]E[(i g2+ p]cov(a 7 j + p, a^ +p ) + cov(d gi+p , d g2+ p)E[a 7 ^ + p]E[a 7 ^ + p] 



Therefore, using (|3. 14[) . (13. 2[) and (|2.1ip . and noting that hi^'i) = h(li) and similarly lj(7 2 ) = ^(72) we get (cf. 
Remark [5]), 



p 1 ™ pa (k 1+ k 2 -2) C0¥ ( fl 7i+P. «72+p) 



CT d ]lfc<o m i h (7i) TO i,(72) = CT d TO 2 +i 2 if 91 = .92 
otherwise. 



Return with the results of Cases 1, 2 and 3 to (|3.10p . In computing the limit of (|3.8p . realize that we need to 
worry about only the case j even and partitions A of pairs. Then T(k : A) is one-to-one with x {ij}£A^(ki, kj). 
Returning to (|3.7p . a moment of thinking gives that 

E[(S n (P)Y] =|pairsof{l,2,...,j}| ( £ w kWl D{k,l) 

\k,l>i 

Since |pairs of {1, 2, . . . ,j}\ = 2 ]/^(j/2) \ ' which are the even moments of the normal N(0, 1), the rest follows. □ 

Corollary 2. Assume that for < e < a b n /n a = 1 + Z n /n e , where lim^^co Z n = Z is a random variable with 
finite moments, mean and variance a\. Then C(k,l) — kla\ and 

D ( k > = \ ^FFIJ ((k-i)/2) ((i-i)/2) J / £ = a and fc > i odd 
[0 otherwise. 

In particular CLT holds for the model (|1.2p . 
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Proof. To compute C(k,l), just notice that (in moments) 

{b n /n a ) k ~ 1 + kZ n /n € + 0{l/n 2e ), 

from which the formula of C(k, I). 

For the rest, there is only one thing we need to do, namely compute 

E EmtMi*)- 

(7i,72)er(fc,z) i<o 

To carry this out, we fist realize the pairs (71,72) G r(fe, i) by fixing 71 in Tfc and then "sliding" up and down 
another fixed ( G Ti, to justify that 

E E wEE^c+s)- 

iez 7 ier fc cer, qez 

For i and C G T;, 2~Z 9 gz ^»(C + l) — ' an d since there are ( ; , 2 ) paths in T;, the rest follows. □ 

4 A Flavor of Free Probability Theory 

Given independent tridiagonal matrices A-y^ n , A2. n , A^ n , one can ask about the joint distribution in the 
moment sense. More precisely, is it true that (here X u>n — —k^A u<n ) 

lim tr n (X lun X i2 , n . . . X lk , n ) exists for any i\, 12, . . . , i m • 

n — >oo 

The answer is yes, but before we do that we need to give a definition. We think about the set C = {1, 2, . . . , r} 
as a set of colors. Then, for a string of colors c = 12, ■ ■ . , ik), from the set C we define 

Vj. = {7 G Tk : each edge j u , j u +i is colored with i u }. 

For a given color u G C, and a colored path 7 G r|, we define 1^(7) as the number of crossings of the line i + 1/2 
with the steps of 7 colored with u. Similarly IF" (7) is the number of flat steps at level i colored with color u. 



Theorem 4. Assume that for each u G C, the entries of the matrix A u>n satisfy, for some a u > 0, 

ra„ k f or an U k > 



lim 1 

n — >oo 



(b u ,n/n a ") k 



with m Ui o — 1 and 

< 00 for any k > 0. 



supE 

n>l 



I d u n I 



Under these assumptions, if all the entries of the matrices are independent of one another and c = Z2, . . . , ik), 
then 

lim tr n pQ 1>n JQ 2in . ..X ik>n ) = ■ ■ ■ — V TT TT m «-i?(7) 

n->oo ai + Ot2 + ■ ■ ■ + Oik + 1 11 11 * 

7Grj«eCiez 

where the limit is in expectation and also almost surely. 

One possible interpretation of this in term of noncommutative probability theory is the following. Assume 
that (X,<j>), (y,i/j) are noncommutative probability spaces, i.e. X, y are unital algebras over the complex 
numbers and ip : X — > C, ip : y — > C are two linear functionals with <p(l) — "0(1) = 1- Assume o' 1; a' 2 , ■ ■ ■ a[ are 
noncommutative random variables on y such that ifj((a' u ) k ) = m U: k- Then the joint distribution of 01, a%, ... a/ 
is described by 

Ha n a l2 ...a lk ) = { ^■"^ E ^ r * ^ 11 k ^ (*) 

1^0 A: odd. 

Note here that this dependence involves m U! fc for fc odd as well, as opposed to the defining relationship from 
(|2.11| which involved only mk for k even. 
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Take for example the case of just two such random matrices. Rescale things out to have a nicer appearance 
to (f> = (ai + Qi2 + 1)0. Let's take two random variables a, b. Then, <f> of a product of odd length in a, b is 0, 
while for products of even length we have (ignore here the presence of ot\ and ai or rescale the functional <j>), 

^b 2 ) = 2^((6') 2 ) 
4>{ab) = 24>(a')4>(b') 

0(a 4 ) = 2^((a') 2 )+4V>((a') 2 ) 2 

0(6 4 ) = 2^(&') 2 )+40((6') 2 ) 2 

4>{abab) = 2^((a') 2 )iP((b') 2 ) + ^(a'f^{b'f 

<f>(a 2 b 2 ) = 4^((a') 2 )V((&') 2 ) + 2^(a') W) 2 

0(a 3 6) = 2^((a') 3 W) + 4V>((a') 2 )V(<W)- 

From this it's quite clear that, with respect to 0, the moments of a and b alone do not determine their joint 
moments. However, imposing the condition that ip(a 2k+1 ) = ip(b 2k+1 ) = for any k > 0, one can do this. For 
example in this case we have 

4>(ab) = 

<j>(abab) = ^cj)(a 2 )<j>(b 2 ) 

4>(a 2 b 2 ) = 4>{a 2 )<t){b 2 ) 
c/)(a 3 b) = 0. 



Another view at these things is the following. Assume that a' and b' are independent random variables 
and the functional ip is just the expectation. Then we consider sequences of iid random variables {X^jigz and 
{Yi}i£z whose distributions are given by the distributions of a' and b' . Consider then the operators A and B 
given in (|2.19p . Now if X is the algebra of infinite dimensional random matrices like A and the functional (j) 
on the algebra generated by A and B is given by K[P(A,B)o o], for any noncommutative polynomial P in two 
variables, then the joint distribution of A and B is given by (*). 

Returning to the general situation from (*), we want to point out that in the case that the variables a' u 
are symmetric, then the noncommutative joint moments of a±, a%, . . . , a; are given in terms of the individual 
moments of oi, 02, ■ ■ ■ , oj. This follows from the fact that all the m^'s involved in the joint moments have k even 
and according to Remark [2] these can be expressed back in terms of the moments of the variables c^. 

We can call these variable "independent" in a certain way and interpret this fact via the relationship between 
the joint moments and the individual moments of each variable. This can be seen by introducing some kind 
of cumulant and express this "independence" as a property of the cumulant. In the classical or free cases of 
independence^ this corresponds to the simple fact that the joint cumulants are the sum of the cumulants of the 
individual variables. 

Finally if all the moments = 1, then the matrices (Xi iT[ , X2, n , . . . , X r>n ) converges in distribution to 
(S, S, . . . ,S) where S is a semicircular random variable, something not very interesting though but due to the 
fact that the coloring does not play any role here. However if the moments mk are not constant equal to 1, then 
the coloring does play an essential role. 

There is also a fluctuation result in this context as follows. 

Theorem 5. Assume that in addition to the properties in the above theorem we have that for some < e u < ot u , 

lim n 2e "cov((6 u ,„/n a ") fc ,(6 u ,„/n a ")') = C u (k, I), 



\n e ((b u>n /n a ") k -E[(b u<n /n^) k }) 

and 

lim vax(d u , n ) = a - . 2 . 



< 00, V k, m > 1. 



Now, take e — min M= i ! ... )T .(e M ) and set 

S n (ii, 22, ■ ■ ■ j ik) = n e+1 / 2 (tr n (Xi 1}n Xi !2tn . . . Xi ktn ) — E[tr ra (Xi 1)n .Xj 3in . . . Xi k ^ n )]) . 
Then the family {S n (ii, i%, ■ ■ ■ ,ik) '■ 1 5= H, ■ ■ ■ > ik < r } converges to a Gaussian family. 
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5 Remarks and Extensions 



5.1 Still Tridiagonal Models 

There are various ways of extending Theorem [1] and Thcorem[3J We refrained to give it in full because the proofs 
would have been overloaded with unnecessary notations and minor differences. 

In both theorems mentioned here we can replace the sequence n a by any sequence a n which satisfies the 
growth rate condition lim,,^,^ n( Q a " — 1) = a. 

The second extension comes from allowing growth in the diagonal part. Namely if we replace the condition 
(|2.8[) by the condition 

lim E[(d n /n> 3 ) k } = m' k 

n — >oo 

then, if < a, the same conclusion holds in Theorem [1] and the same conclusion under Theorem [3J with 
the condition (|3.2p replaced by lim^oo var(d„/n /3 ) = <7<j. However if [3 = a, then the conclusions still hold, 
nonetheless the formulae of L k become 

Lk= ^F+\ ^ n m ii(7)<(7)> 

7 er fc ur- *<o 

while in Theorem [31 condition (|3.2p has to be replaced by 

lim n 2t cov((d n /n a ) k ,(d n /n a ) 1 ) = C'(k,l) 

n — >oo 

the only difference here is that the convariance matrix D e (k, I) now depends also on C'(k, I) and m' k . 

If /3 > a, the scaling of the matrix X n has to be changed to X n — -^A n . The conclusions of both theorems 
hold with the appropriate changes since now the dominating terms are the ones on the diagonal. For example 
the (|2.1ip . becomes 

Lk = WTi m ' k - 

We leave to the reader to see how the changes in Theorem [3J have to be done. 

Another extension is obtained by dropping the independence of the entries. We can replace that in Theorem[T] 
by a more relaxed version. 

Remark 6. Assume that for any 7 G T k U T k , there is a number m 7 so that for a certain a > 0, 

lim -T^ E l a l+n] = m ~f- 

rwoo n 



Then if X n = -^A n 

lim tr n (X^) = L k 

n — >oo 

where L k is computed by 

L k = — >^ m~. 

7er fc ur fc 

In particular one can apply this to the cases when the matrix A is obtained from another tridiagonal matrix B n 
which has independent entries by replacing each entry with a function of the other nearby entries in a finite 
range. For example one can replace the nonzero entries in B n by the average of the neighbors nearby it in a 
finite range. Another example is the Laguerre (3 models discussed in f3j] and J^, or more general the models in 
which each entry in B is replaced by a polynomial in the variables lying in finite neighborhood of the entry. 



5.2 Band Diagonal Models 

We can extend the results so far to a more general setting by allowing not only one subdiagonal but more than 
one. In this case we take symmetric matrices of the form A n = so that ai.j = for \i — j\ > w, where 

w is the width of the band and all entries are independent. Denote b v ,i — ai^ v +i. 

In this case we can consider the problem of convergence of the eigenvalues and of the fluctuations. Before we 
give this extension, let us define the needed objects. 
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Set 

Tk, w = {7 = (*i>*2, ■ ■ ■ ,h+i) € Z fc+1 : i'l = «fc+i, |»« " < w,max(7) = 0}. 

Then we define for any path A, luj] (A) to be the number of steps i u , i u +\ so that {i u , i u +i} = {i, j}- In particular, 
for the notations we already used we get lj(A) = (A) and ffi(A) = l[ i] ,](A). 

Notice here the equivalent of the formula (|2.6|) as 

ax = n ^r i(A) - 

ieZ,0<v<w 

Theorem 6. Assume that for each < v < w, and given a v > 0, there is m„,fc so that 

|lim, woo E[(a ij /n Q ") fc ] = m V:k , if\i~j\ = v,a v = a 
[sup„E[(a ij /n Q ")' c ] < 00, if\i~j\ = v,a v < a 

with < a = max(a„ : 1 < v < w) and the convention that m v fi = 1 for any v. Then, for X n = ^A n , one has 
that 

lim tr n (Z*) = L fc 



both in average and almost surely. Moreover, 



L k = 



k= 7^TT s n m i[M + ,](7). 

76T fc , ra i<0,0<t)<tu 

where 

rrik if ol v = a 



if Olv < Q!- 



This theorem says that in fact those subdiagonals not scaled by the maximum power n a , do not contribute 
to the limit Lfc- 

Let us point out that one can extend this to a statement in which the independence condition is dropped 
and one gets a version of Remark [5J 

Similar versions of the first part of Proposition [1] can be proved in this context too. Namely, if each of the 
moments (mv,k)kLi come from the moments of a compactly supported measure, then the moments Lk also come 
from a compactly supported measure. In addition, if there are some numbers m v so that m v _k = rn v for all 
k > 1, then the corresponding measure with the moments given by Lk is the distribution of J2™=i m v{z v + z~ v ) 
under the Haar measure of the circle S . 

There is also a version of Theorem [3l 

Theorem 7. In addition to the conditions given in the above Theorem, assume that for each < v < w, there 
is < e v < a v so that for any k, I > 1 

lim n 2£ "cov((6 !J ,„/n^) fc ,(6 !J ,„/n Q ") fc ) = C u (k,l) 

n — >oo 

and 



E 



< 00, V k, m > 1, 



\n^((b Vin /n^) k -E[(b v , n /n a ") k })\ 
Let e = min(e„ : < v < w) and define 

S n (fc)=n e + 1 /2( t r n (^)-E[tr n (^)]). 
Then the family {S n (k)}' k x L 1 converges to a family of Gaussian random variables. 



References 

[1] Greg W. Anderson and Ofer Zeitouni, A CLT for a Band Matrix Model, Probability Theory and Related 
Fields 134 (2006), no. 2, 283-338. 

[2] Z.D. BAI, Methodologies in Spectral Analysis of Large Dimensional Random Matrices: a review, Statistica 
Sinica 9 (1999), no. 3, 611-661. 



25 



[3] Ioana Dumitriu and Alan Edclman, Matrix Models for Beta Ensembles, Journal of Mathematical Physics 
43 (2002), no. 11, 5830-5847. 

[4] , Global Spectrum Fluctuations for the (3-Hermite and (3-Laguerre ensembles via matrix models, 

Journal of Mathematical Physics 47 (2006), no. 11, 5830-5847. 

[5] P Flajolet, Combinatorial Aspects of Continued Fractions, Discrete Mathematics (1980), no. 32, 125-161. 

[6] K. Johansson, On Fluctuations of Random Hermitian Matrices, Duke Math. J. 91 (1998), 1-24. 

[7] Madan Lai Mehta, Random Matrices, second edition ed., Academic Press, 1991. 

[8] James A. Mingo and Roland Spcicher, Second Order Freeness and Fluctuations of Random Matrices: I. 
Gaussian and Wishart matrices and Cyclic Fock spaces, preprint math. OA/0405191 (2005). 

[9] A. Soshnikov, Universality at the Edge of the Spectrum in Wigner Random Matrices, Comm. Math. Phys. 
207 (1999), no. 3, 697-733. 

[10] A. Soshnikov and Y Sinai, Central Limit Theorem for Traces of Large Symmetric Matrices with Independent 
Matrix Elements, Bol. Soc. Bras. Mat. 29 (1998), 1-24. 

[11] Roland Spcicher, Combinatorial Aspects of Free Probability Theory, Lecture notes, 2005, "Freie Wahrschcin- 
lichkeitstheorie" , Goettingen. 

[12] Daniel W. Stroock, Probability Theory, an Analytic View, revised ed., Cambridge University Press, 2000. 

[13] Gerald Teschl, Jacobi Operators and Completely Integrable Nonlinear Lattices, AMS, 1999. 

[14] Dan Voiculescu, Limit Laws for Random Matrices and Free Products, Invent, math. 104 (1991), 201-220. 

[15] Voiculescu, D., Dykcma, KJ, and Nica, A., Free Random Variables, "CRM Monograph Scries, vol. 1, AMS, 
1992. 

[16] E. P. Wigner, Characteristic Vectors of Bordered Matrices with Infinite Dimensions, Ann. Math. 62 (1955), 
548-564. 

[17] , On the Distribution of the Roots of Certain Symmetric Matrices, Ann. Math. 67 (1958), 325-327. 



26 



